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Abstract 

We investigate exclusive nonleptonic B meson decays B -^ (j)K in perturba- 
tive QCD formalism. It is shown that the end-point (logarithmic and linear) 
singularities in decay amplitudes do not exist, after kx and threshold resum- 
mations are included. Power counting for emission and annihilation topologies 
of diagrams, including both factorizable and nonfactorizable ones, is discussed 
with Sudakov effects taken into account. Our predictions for the branching 
ratios B{B ^ 0^) ~ 10 x 10^^ are larger than those (~ 4 x 10^^) from the 
factorization approach because of dynamical enhancement of penguin contri- 
butions. Whether this enhancement is essential for penguin-dominated modes 
can be justified by experimental data. 



PACS number(s): 12.38.Bx, 13.35. Hw, 12.38.Qk, ll.lO.Hi 



I. INTRODUCTION 

Perturbative QCD (PQCD) factorization theorem for the semileptonic decay B -^ tiW 
has been proved in |jl]], which states that the soft divergences in the B ^ n form factor can 
be factorized into a hght-cone B meson distribution amphtude and the colhnear divergences 
can be factorized into a pion distribution amphtude order by order. The remaining finite 
contribution is assigned into a hard amphtude, which is calculable in perturbation theory. 
A meson distribution amplitude, though not calculable, is universal, since it absorbs long- 
distance dynamics, which is insensitive to specific decays of the b quark into light quarks 
with large energy release. The universality of nonperturbative distribution amplitudes is one 
of the important ingredients of PQCD factorization theorem. Because of this universality, 
one extracts distribution amplitudes from experimental data, and then employ them to 
make model-independent predictions for other processes. In this paper we shall assume that 
PQCD factorization theorem holds for two-body nonleptonic B meson decays, to which there 
is no difficulty to generalize the proof in |l[. The one- loop proof for the PQCD factorization 
of two-body decays has been given in 0. 

The PQCD formalism for the charmed decays B -^ D^*^7r(p) [Q,^ is restricted to twist-2 
(leading-twist) distribution amplitudes. For charmless decays such as B ^ Kn, tttt and KK 
IQ-^, contributions from two-parton twist-3 (next-to-leading-twist) distribution amplitudes 
are introduced via the penguin operators O^s in the effective Hamiltonian for weak decays. 
It has been argued [^ that two-parton twist-3 contributions are in fact not suppressed 
by a power of 1/Mb, Mb being the B meson mass. A compelte leading-power PQCD 
analysis of the heavy-to-light B ^ n, p form factors, including both twist-2 and twist-3 



contributions, has been performed in ^^. There exist many other higher- twist sources in 
B meson decays, whose contributions are indeed down by a power of I/Mb- These sources 
include the B meson and b quark mass difference A = Mb — rrib, the light quark masses 
mu, TTLd and rus, and the light pseudo-scalar meson masses M^r and Mk- Those from three- 
parton distribution amplitudes are further suppressed by the coupling constant a^. All these 



sub-leading contributions will be neglected in the current formalism. 

In this work we shall perform a PQCD analysis of the B -^ (j)K decays up to corrections 
of 0{A/Mb)- These modes involve different topologies of diagrams, such as factorizable 
(also nonfactorizable) emission and annihilaton. We predict the branching ratios and CP 
asymmetries of the modes 



B', 



(j)K' 



B' 



(j)K^ . 



(1) 



It will be found that two-parton twist-3 contributions are comparable to twist-2 ones as 
expected. Our predictions for the braching ratios 



10.211?) X 10"^ 



B{B^ -^ <pK^) 



(2) 



are larger than those from the factorization approach [11] and from the QCD factorization 



approach |12|, which are located within 4.3tf;4 x 10 ^ [13,n] with the uncertainty arising 



mainly from the inclusion of annihilation contributions [Q. The mechanism responsible for 
the larger branching ratios in the PQCD formalism is dynamical enhancement of penguin 
contributions MM- Note that the current experimental data of B{B^ — > (pK^)., 



CLEO [15] : 



Belle 16 



BaBar 17 : 



(5.5lf:^±0.6) X 10" 



10.6lf:^±2.2) X 10" 
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and those of 5(5°^ 0ir° 



CLEO [15] : 
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BaBar [IT? 



< 12.3 X 10" 
(8.7l|^ ± 1.5) X 10"^ 
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are still not very consistent with each other. 



The PQCD factorization formulas for B -^ (pK decay amplitudes have been derived 
independently in |T^. Here we shall further discuss the power behavior of the factorizable 
emission and annihilation amplitudes, and the nonfactorizable amplitudes in 1/Mb, whose 
relative importance is given by 

emission : annihilation : nonfactorizable = 1 : : , (5) 

Mb Mb' ^ ^ 

with mo being the chiral symmetry breaking scale. In the heavy quark limit the annihilation 
and nonfactorizable amplitudes are indeed power-suppressed compared to the factorizable 
emission ones. Therefore, the PQCD formalism for two-body charmless nonleptonic B meson 
decays coincides with the factorization approach as Mb — > C)0. We shall also explain why 
the annihilation amplitudes are mainly imaginary and investigate theoretical uncertainty in 
the PQCD approach. 

Dynamical enhancement is a unique feature of the PQCD approach, which does not exist 
in the factorization or QCD factorization approach. We argue that the B —>■ (j)K modes are 
more appropriate for testing this mechanism in penguin-dominated nonleptonic B meson 
decays compared to the B -^ Kit decays |^. The large B — > Kit branching ratios may 
not be regarded as an evidence of dynamical enhancement: they can also be achieved by 
chiral enhancement (the kaon is a pseudo-scalar meson) and by choosing a large unitarity 
angle 03 ~ 120° [0, which leads to constructive interference between penguin and emission 
contributions. The B -^ (pK modes are not chirally enhanced, because is a vector meson, 
and insensitive to the variation of the angle 03, because they are pure penguin processes. If 
the data of B{B -^ (j)K) are settled down at values around 10 x 10~^ in the future, dynamical 
enhancement will gain a strong support. 

On the other hand, precise measurement of the CP asymmetry in the B —>■ (j)K decays 
is important for new physics search and for the determination of the unitarity angle 0i 
with high degree of accuracy pO|j2T[] . This measurement is experimentally accessible at 



the early stage of the asymmetric B factories. The B -^ (pK decays arise from penguin 
(loop) effects, while the B -^ J /ipK decays involve only tree amplitudes. The search for 



different CP asymmetries in the B -^ J/ipKs and (pKs decays, with the common source 
from B^-B^ mixing, provides a promising way to discover new physics p2| , p3| : a difference 
of \Acp{J /ipKs) — Acp{(t>Ks)\ > 5% would be an indication of new physics. This subject 
will be addressed elsewhere [Q. Furthermore, the and J/ip mesons are bound states of 



the ss and cc quarks. It is also interesting to compare their branching ratios, which reflect 
the mass effect from charm quarks. 

We demonstrate the importance of kr and threshold resummations by studying the B — >■ 
TT, K transition form factors in Sec. II. The power counting and the factorization formulas for 
various topologies of amplitudes are given in Sec. III. The numerical analysis is performed in 
Sec. IV. Section V is the conclusion. Twist-2 and two-parton twist-3 distribution amplitudes 
for the kaon and for the meson are defined in the Appendix. 

II. SUDAKOV SUPPRESSION 

In this section we briefly review the importance of k^ and threshold resummations for an 
infrared finite PQCD calculation of heavy-to-light transition form factors |TD|. Consider the 



leading diagram shown in Fig. 1 for the B —>■ (pK decays in the kinematic region with a fast- 
recoil kaon. The B meson momentum Pi, the meson momentum P2 and the longitudinal 
polarization vector e, and the kaon momentum P3 are chosen, in light-cone coordinates, as 

^ Mb, . Mb, n , Mb, 9 . 

Pi = -^(1,1,0t), P2 = ^(l,r5,OT), P3 = ^(0,l-r5,OT), 

.^^(1,-J,0.), (6) 

with the ratio r^ = M^/Mb, M^ being the meson mass. At the end of the derivation 
of the factorization formulas, the terms r^ ~ 0.04 in the above kinematic variables will be 
neglected. We treat the kaon as a massless particle, and define the ratio tk = mo/MB for 
the kaon, which will appear in the normalization of the twist-3 kaon distribution amplitudes. 
The B meson is at rest under the above parametrization of momenta. 



0(A,e) 




kx d, u ^3 

FIG. 1: Leading contribution to the B —* (j)K decays, where H denotes the hard 
amphtude and ki, i = 1, 2, and 3, are the parton momenta. 

A. kx and threshold resummations 

It has been known that the lowest-order diagram in Fig. 2(a) for the B -^ tt form factor 
F^"^ is proportional to l/(xiX3) without including parton transverse momenta kx, where 
Xi = ki / Pi [x-i = k^ / P^) is the momentum fraction associated with the spectator quark 
on the B meson (pion) side. If the pion distribution amplitude vanishes like x^ as X3 — > 
(in the twist-2 case), F^'^ is logarithmically divergent. If the pion distribution amplitude 
is a constant as X3 — ;> (in the twist-3 case), F^'^ even becomes linearly divergent. These 
end-point singularities have caused critiques on the perturbative evaluation of the 5 — > vr 
form factor. Several methods have been proposed to regulate the above singularities. An 



on-shell h quark propagator has been subtracted from the hard amplitude as X3 — » in p5 
However, this subtraction renders the lepton energy spectrum of the semileptonic decay 
B — > -kW vanishes as the lepton energy is equal to half of its maximal value. Obviously, 
this vanishing is unphysical, indicating that the subtraction may not be an appropriate way 
to regulate the singularity. The subtraction also leads to a value of F^"^ at maximal recoil, 
which is much smaller than the expected one 0.3. A lower bound of 0:3 of ©(A/M^) has been 
introduced in [^ to make the convolution integral finite. However, the outcomes depend 
on the cutoff sensitively, and PQCD loses its predictive power. 
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(a) 
FIG. 2: Lowest-order diagrams for the B —>■ tTjK transition form factors. 




A self-consistent prescription has been proposed in ^^ , where parton transverse moemta 
kT are retained in internal particle propagators involved in a hard amplitude. In the end- 
point region the invariant mass of the exchanged gluon is only O(A^) without including 
kx- The inclusion of kj- brings in large double logarithms a^ In (kr/MB) through radiative 
corrections, which should be resummed in order to improve perturbative expansion, kx 
resummation |28,29| then gives a distribution of kj- with the average (k^) ~ 0{AMb) for 
Mb ~ 5 GeV. The off-shellness of internal particles then remains of 0{AMb) even at the 
end point, and the singularities are removed. Hence, it can not be self-consistent to treat 
kx as a higher-twist effect as the end-point region is important. The expansion parameter 
as{AMB)/'n' ~ 0.13 is also small enough to justify PQCD evaluation of heavy-to-light form 
factors 0,0. This result is so-called Sudakov suppression on the end-point singularities in 
exclusive processes PU| . 

The above discussion applies to the B ^ K form factor and to the B —>■ (pK decays, 
fcy resummation of large logarithmic corrections to the B, cf) and K meson distribution 
amplitudes lead to the exponentials Sb, Sa, and Sk, respectively: 



Ssit) = exp 
S^{t) = exp 
Sxit) = exp 



Ji/bi n 



Jl/b2 U 
-S{X3P^, 63) - s((l - X3)P^, 63 



2/' ^liasm')) 

l/bs fi 



(7) 



with the quark anomalous dimension 7 = — a^/vr. The variables 61, 62, and 63, conjugate 
to the parton transverse momenta kij-, k2T, and k^^, represent the transverse extents of the 
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5, and K mesons, respectively. The expression for the exponent s is referred to p8|-|30|. 
The above Sudakov exponentials decrease fast in the large h region, such that the B —>■ (pK 
hard amplitudes remain sufficiently perturbative in the end-point region. 

Recently, the importance of threshold resummation [pl| -|33[] has been observed in exclu- 
sive B meson decays EM. As X3 ^ (to be precise, X3 ~ 0{K/Mb)) in Fig. 2(a), the internal 



h quark, carrying the momentum Pi — k^, becomes almost on-shell, indicating that the end- 
point singularity is associated with the b quark. Additional soft divergences then appear at 
higher orders, and the double logarithm a^ In X3 is produced from the loop correction to 
the weak decay vertex. This double logarithm can be factored out of the hard amplitude 
systematically, and its resummation introduces a Sudakov factor Sti^x^) into PQCD factor- 
ization formulas ^^. Similarly, another lowest-order diagram Fig. 2(b) with a hard gluon 



exchange between the d{s) quark and the spectator quark gives an amplitude proportional to 
l/(xfx3). In the threshold region with xi -^ (to be precise, xi ~ 0(A^/M^)), additional 
collinear divergences are associated with the internal d{s) quark carrying the momentum 
P3 — fci. The double logarithm a^ln^xi is then produced from the loop correction to the 
weak decay vertex. Resummation of this type of double logarithms leads to the Sudakov 
factor St{xi). 

The above formalism can be generalized to factorizable annihilation diagrams easily. For 
the lowest-order diagram with the internal quark carrying the momentum P2 + k^, the end- 
point region corresponds to X3 — ^ 0. Hence, threshold resummation of the double logarithm 
gives the Sudakov factor Stix^). For the lowest-order diagram with the internal quark 
carrying the momentum P^ + k2, the end-point region corresponds to X2 — > 0. Threshold 
resummation of the double logarithm then gives the Sudakov factor St{x2). The Sudakov 
factor from threshold resummation is universal, independent of flavors of internal quarks, 
twists and topologies of hard amplitudes, and decay modes. To simplify the analysis, we 
have proposed the parametrization pO|| , 

S,{x) = ^-^^^^^l^[x{l-x)r. (8) 

VttI (1 + c) 



with the parameter c = 0.3. This parametrization, symmetric under the interchange of x 
and 1 — x, is convenient for evaluation of factorizable annihilation amplitudes. It is obvi- 
ous that threshold resummation modifies the end-point behavior of the meson distribution 
amplitudes, rendering them vanish faster at a; ^ 0. 

Threshold resummation for nonfactorizable diagrams is weaker and negligible. The rea- 
son is understood as follows. Consider the diagram with a hard gluon exchange between the 
spectator quark and the s quark in the meson, in which the internal s quark carries the 
momentum k2 — ki + k^. To obtain additional infrared divergences, i.e., the double loga- 
rithms, the s quark must be close to mass shell. We then have the threshold region defined, 
for example, by ki ~ 0{A'^/Mb), k2 ~ 0{h? /Mb)., and k-^ ~ 0{Mb) simultaneously. That 
is, this region has more limited phase space compared to that for factorizable amplitudes. 
Furthermore, soft contribution to a pair of nonfactorizable diagrams cancels, such that the 
end-point region is not important. Based on the above observations, we shall not include 
threshold resummation for nonfactorizable amplitudes. 

kT and threshold resummations arise from different subprocesses in PQCD factorization. 
They can be derived in perturbation theory, and are not free parameters. Their combined 
effect suppresses the end-point contributions, making PQCD evaluation of exclusive B me- 
son decays reliable. If excluding the resummation effects, the PQCD predictions for the 
B ^ K form factor are infrared divergent. If including only kx resummation, the PQCD 
predictions are finite. However, the two-parton twist-3 contributions are still huge, so that 
the B ^ K form factor has an unreasonably large value F^^ ~ 0.57 at maximal recoil. 
The reason is that the double logarithms a^ In^ x have not yet been organized. If including 
both resummations, we obtain the reasonable result F^^ ~ 0.35. This study indicates the 
importance of resummations in PQCD analyses of B meson decays. In conclusion, if the 
PQCD evaluation of the heavy-to-light form factors is performed self-consistently, there exist 
no end-point singularities, and both twist-2 and twist-3 contributions are well-behaved. 

The mechanism of Sudakov suppression can be easily understood by regarding a meson 
as a color dipole. In the region with vanishing kx and x, the meson possesses a huge extent 
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in the transverse and longitudinal directions, respectively. That is, the meson carries a 
large color dipole. At fast recoil, this large color dipole, strongly scattered during B meson 
decays, tends to emit infinitely many real gluons. However, these emissions are forbidden 
in an exclusive process with final-state particles specified. As a consequence, contributions 
to the B ^ 7T,K form factors from the kinematic region with vanishing k^ and x must be 
highly suppressed. 

B. Form Factors 

We demonstrate that kx and threshold resummations must be taken into account in 
order to obtain reasonable results for the B ^ K,tt form factors. In the PQCD approach 
these form factors are derived from the diagrams with one hard gluon exchange as shown in 
Fig. 2. Soft contribution from the diagram without any hard gluon is Sudakov suppressed 
||1CI|| . For a rigorous justification of this statement in QCD sum rules, refer to [^. The two 
form factors F^Q^q^) are defined by 



{K{p,)\m^,smB{Pi)) 



■^BK/2\ 



(^1 + ^3). - \,^ % 



+Fi^{<l')^^^^% , (9) 



where g = Pi — P3 is the outgoing lepton-pair momentum. F^^ and F^^ at maximal recoil 
are, quoted from Eq. ( |39D below, written as 

F^^{0) = 87tCfMI dxidxs 6irf6iM&3$B (a;i, &i) 
Jo Jo 

X I [(1 + X3) $X (X3) + (1 - 2X3) TK ($^ (X3) + $J, (X3))] 
Xa,(t«)5B(ti'))5;,(t«)/le {xuxs, 61, 63) 

+2rx$^ (X3) as{tf^)SB{tf'>)SK{tf^)he {x^, x,, h, b,) | , (10) 

with Cp = 4/3 being a color factor. The hard function he{xi, x^, 61, 63), referred to Eq. 
contains the threshold resummation factor Stix^). The hard scales tg are defined in Eq. 
The expression for the B ^ tt form factor F^'^{q^) is similar to Eq. (p!oD . 
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In Eq. (1^) we have included the complete two-parton twist-3 distribution amplitudes 
associated with the pseudo-scalar and pseudo-tensor structures of the kaon. We adopt 
the mass mo = 1.7 GeV @J§], the B meson distribution amplitude $b proposed in |^,|^ 
(see Eq. (0)), and the kaon distribution amplitudes ^k, '^^, and $^ derived from QCD 
sum rules [^ (see Eqs. (p5|)-(p7|)). Since $b is still not well determined, we consider 
the variation of its shape parameter within 0.36 GeV < ujb < 0.44 GeV ^^ around its 
central value ujb = 0.4 GeV. Turning off kx and threshold resummations and fixing ««, 
the B ^ 7i,K form factors are divergent and not calculable. With fc^ resummation, the 
results become finite as shown in Table I, but still much larger than the expected ones 0.3- 
0.4. Further including the threshold resummation effect, we obtain the reasonable values 
F^^(O) ~ 0.35 ±0.06 and F^^(O) = 0.30 ±0.04. These ranges of the form factors have been 
usually adopted as model inputs in the literature, and are consistent with the results from 
lattice calculations p8|- ^] extrapolated to the small g^ region and from light-cone QCD 



sum rules ET|,H2 



We also present our results of the time-like form factors F"^^, which govern the factor- 
izable annihilation amplitudes. The factorization formulas are, quoted from Eqs. (^Tj) and 
(0), given by 

^(v-A) = 8itCfMI / dx2dx3 / h2dh2hdh'i 

X I [(1 - Xg) ^4,{X2)^K {x-i) + 2rKr^% (X2) ((2 - X3)$^ (Xg) ± X3$^ (xg)) 
Xa,(t«)5^(t«)5^(t«)/i, {X2, 1 - X3, 62, h) 

- X2<^^ {X2) <^K (xg) ± 2rxr^ ((1 ± X2)% (xa) - (1 - Xs)^^ (xs)) $^ (xg) 

Xa,(ti2))5^(t(2))5^(t(2))/,^ (1 _ X,, X2, 63, &2) I . (11) 



(V+A) 



/•I roo 

—SnCpM'^ / dx2dx3 / b2db2bsdbs 
Jo Jo 

X I [2ri^(l - X3)$0(X2) ($^(3^3) ± $x(x3)) ± ^r^% (X2) ^k (X3) 
xa,(t«)5^(t«)5x(t«)/i, (X2, 1 - X3, 62, bs) 
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+ 



4ri^$<^(x2)$^(x3) + 2x2r^ (^^(xs) - $*^ (xs)) ^k{.x^) 



xa.M'^)SM'^)SK{t'^)ha (1 - xa, X2, 63, h) | . (12) 

The hard function ha{x2, X3, 62, ^3), referred to Eq. (^), contains the threshold resummation 
factor Sti^xs). The hard scales ta and the meson distribution amplitudes ^^, $0 and ^^ 
are defined in Eq. ( ^7\ ) and in Eqs. (|5^-(|51D, respectively. It is obvious that Ffy_^-. con- 
tains both logarithmic and linear divergences, and Ffy_^^-. contains logarithmic divergences 
|l^,43|, if the Sudakov factors are excluded. To include annihilation contributions in the 



QCD factorization approach, several arbitrary complex infrared cutoffs must be introduced 



to regulate the above end-point singularities |Q. These cutoffs are process-dependent. 
Hence, the PQCD approach with the Sudakov effects has a better control on annihilation 
contributions. 

We obtain Ffy_^-^ = (1.78 + 0.63i) x lO'^ ^^^ F({,^^) = (-3.48 + 13.52i) x lO^^ for 
cub = 0.4 GeV and mo = 1.7 GeV. The smaller value of Ffy_^'. is due to mechanism similar 
to helicity suppression: note the partial cancellation between the two terms in Eq. {^A\j. 
-^(y+yi) is mainly imaginary, whose reason will become clear after we explain the power 
behavior of the various decay amplitudes in the next section. The time-like form factors are 
independent of the shape parameter ub- The change of the B —>■ tt,K form factors and of 
the time-like form factors with the chiral enhancement factor mo is displayed in Table II. 

III. POWER COUNTING AND DECAY AMPLITUDES 

We discuss the power counting rules in the presence of the Sudakov effects and present 
the factorization formulas for the B —>■ (pK decays. The effective Hamiltonian for the flavor- 
changing b -^ s transition is given by 



Grp 

V ^ q=u, 



^eff = ^= I] Vg 



q=u,c 
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i=3 



(13) 



with the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vq = V*gVqh and the opera- 
tors 

13 



Oi"^ = {Siqj)v-A{qjbi)v-A , Oi"^ = {Siqi)v-A{qjbj)v-A , 

O5 = {Sibi)v-Aj2i^j1j)v+A , Oe = {sibj)v-A^{qjqi)v+A , 

3 3 

Oi = -isibi)v-A^eg{qjqj)v+A , Os = -{sibj)v-Aj2^g(^j1i)v+A , 

3 3 

O9 = -(Si6i)y_A Yl '=^giQj1j)v-A , Oio = -{Sibj)v-A Yl (^qiQj(li)v^A , 



(14) 



i and j being the color indices. Using the unitarity condition, the CKM matrix elements for 
the penguin operators O3-O10 can also be expressed as V^ + Vc = —Vt- The unitarity angle 
03 is defined via 



Vub = |Kfe|exp(-z03) . 

Here we adopt the Wolfenstein parametrization for the CKM matrix upto 0{X^) 
iVud Vus KA / I-AV2 A A\\p~iri)\ 



Vcd Vcs Vcb 

V Vtd Vts Vtb ) 



-A 1 - AV2 

\A\^{\-p-i^) -A\^ 



A recent analysis of quark- mixing matrix yields |^| 



A\^ 
1 



(15) 



(16) 



A = 0.2196 ±0.0023, 
A = 0.819 ±0.035, 
i?6 = Jp^ + r/2 = 0.41 ± 0.07 



(17) 



The hard amplitudes contain factorizable diagrams, where hard gluons attach the valence 
quarks in the same meson, and nonfactorizable diagrams, where hard gluons attach the 
valence quarks in different mesons. The annihilation topology is also included, and classified 
into factorizable and nonfactorizable ones according to the above definitions. The B — > (\)K 
decay rates have the expressions 



r= ^ ^ \A? 

327r ' ' 



(18) 
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The amplitudes for Bd -^ 4'K^, Bd — > 'pK'^ ■, B^ -^ ipK^ and B -^ (pK are written, 
respectively, as 

A = UV;Ff^^^ + V^M^'^'^ + fBV;F[^''^ + y/A^r^"^) , (19) 

A = UVtF[^'^ + VtM^^'^ + fsYtF^^'^ + \4-Mr(") , (20) 

A+ = /^\/;Ff w + \/;A<f (^) + /Bi^;Ff («) + v:m^J-'^ - fBv:Fa - v:Ma , (21) 

A- = UVtF^^'^ + ^t^f (^) + /B^iFf (") + V^iA^f (") - fBVuFa - VuMa . (22) 

In the above expressions F {A4) denote factorizable (nonfactorizable) amplitudes, the sub- 
scripts e (a) denote the emission (annihilaiton) diagrams, and the superscripts P{q) denote 
amplitudes from the penguin operators involving the q-q quark pair, and fs {f<f,) is the B 
(0) meson decay constant. 

A. Power Counting 

Before presenting the factorization formulas of the above amplitudes, we discuss their 
power behavior in I/Mb- The spectator quark in the B meson, forming a soft cloud around 
the heavy b quark, carries momentum of 0(A). The spectator quark on the kaon side carries 
momentum of 0{Mb) in order to form the fast-moving kaon with the s quark produced in 
the b quark decay. These dramatic different orders of magnitude in momenta explain why 
a hard gluon is necessary. Based on this argument, the hard gluon is off-shell by order of 
AMb- This special scale, characterizing heavy-to-light decays, is essential for developing 
the PQCD formalism of exclusive B meson decays. Below we shall explicitly show how to 
construct this power behavior, and argue that all the topologies of diagrams should be taken 
into account in the leading-power PQCD analysis. 

We start with the twist-2 contribution to the factorizable emission amplitude, which 
contains the integrand (see Eq. (|TU])), 

[xiX3M| + 0(AMb)][x3M| + 0(AMb)] ' ^ ^ 
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where the factor X3 in the numerator conies from the twist-2 kaon distribution amphtude 
^kI^s), the first and second factors in the denominator come from the virtual gluon and 
quark propagators, respectively. The terms of 0{AMb) represent the order of magnitude of 
parton transverse momenta k^ under Sudakov suppression. The momentum fraction Xi is 
assumed to be of 0{A/Mb). Here we concentrate on the important end-point region with 
X3 ~ 0{A/Mb). It is trivial to find that in the end-point region P^ behaves like 

r^ ~ -^ , (24) 

AMb ' ^ ^ 

as argued above. 

We then consider the first twist-3 term in Eq. (0), 



re3 



tkMI 



[xiXsMi + 0(AMb)][x3M| + O(AMb)] 
For small 0:3 ~ 0{A/Mb), we have the power law, 



(25) 



3 tk mo 1 



A2 A AMb ^ ' 

Hence, the twist-3 contribution is not power-suppressed compared to the twist-2 one in the 
Mb — ^ 00 limit [P|,pi|, though the two-parton twist-3 distribution amplitudes are propor- 
tional to the ratio tk = rrio/MB- The presence of the potential linear divergence in the 
twist-3 contribution modifies the naive power counting rules from twist expansion of meson 
distribution amplitudes. The power behavior of the other twist-3 term in Eq. ([T0| ) is the 
same. 

A folklore for annihilation contributions is that they are suppressed by a power of the 
small ratio Jb/Mb, and negligible compared to emission contributions. The annihilation 
amplitudes from the operators Oi^2,3,4 with the structure {V — A)(y — A) are small because 
of mechanism of helicity suppression. This argument applies exactly to the B ^ titt de- 
cays, and partially to the B -^ Kn and (j)K decays, since the kaon distribution amplitudes 
are not symmetric in the momentum fraction x, and the two final-state mesons are not 
identical. Those from the operators O^^ with the structure {S — P){S + P) survive under 
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helicity suppression @J§]. Below we shall argue that the annihilation amplitudes from O^^e 
are proportional to 2rK, which is in fact of 0(1) for Mb ~ 5 GeV. That is, annihilation 
contributions vanish in the heavy quark limit nib -^ oo, but are important for the physical 
mass nib- 

Referred to Eq. (0), the integrand for the twist-2 annihilation amplitudes from Oi^2,3,4 
is written as, 



ra2 _ X2XIMI 

[X2X^MI - 0{KMb) + ie][x^Ml - 0{KMb) + le] 






where the factor X2 in the numerator comes from $(/,(a;2) and x| comes from x^^k{xz)- Here 
we have interchanged x^ and 1 — x-^ for the convenience of discussion. In the region with 
X3 ~ 0{A/Mb) and with arbitrary X2, we obtain the power laws of the real and imaginary 
parts of Eq. (|3), 

Mli^-A))-^, M/(t-A))-T^- (28) 

B B 

Hence, the twist-2 annihilation contributions from Oi 2,3,4 are negligible. 
The integrand for the Olr"^) terms in Eq. (^TJ) is written as. 



ra4 _ 2rKr^M^B 

[X2X3MI - O(AMb) + «e][x3M| - O(AMb) + te] 



Ta4: _ z,, Ki<t>iy^B /0Q\ 

Hv-A) - TZr~rT72 7T7T7FTT~~TTr"T72 TTTTTTTT""? • V^^J 



We express the quark propagator as 



1 P ZTT / „ / A 



M ^3 - O -— , (30) 



X3MI - O(AMb) + ze X3MI-0{AMb) M| V"'^ \Mb , 

where P denotes the principle-value prescription. The gluon propagator is expressed in a 
similar way. The real part then behaves like 



X2X3MI - 0{KMb)] [x^Ml - 0{AMb) 
(I 7r2 \ Ms I ( n^A\ , , 
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The imaginary part of Eq. (^) behaves hke 

A 



/m(/(V_^)) = -2-KrKr^ 



S[xs-0 



[x2XsMl - 0{AMb)] \ \Mb 






[x^MI-0{KMb)] \ ' ' \Mb) , 

2tX Ma, 271 

For the above estimates the end point 0:3 ~ 0{A/Mb) with arbitrary X2 corresponds to 
the important region. Note that the ratio M^/A is of 0(1). Compared to Eq. (p^), the 
first term in Eq. (0), multiphed by a chiral factor 2rK ~ 0(1), is not down by the small 
ratio Jb/Mb- The second term in Eq. (|31| ) and the imaginary part, though scaling like 
1/M^, are enhanced by the factors vr^ and 27r, respectively. However, the mechanism of 
helicity suppression renders these annihilation contributions turn out to be small (~ 1/M^) 
as shown in Table II. 

Referred to Eq. ([12D , the general integrand for the twist-3 annihilation amplitudes from 
05,6 is given by, 

ra3 ^ 2rKX2X^Ml 

(^+^) [x2X^Ml-0CAMB)+ie]\x^Ml~0CAMB)+it]' ^ ' 

By means of a similar argument, the real and imaginary parts behave, respectively, like 

fle(/,-,,,) ~ 2r.. Jj- (1 - ^) . (34) 



AMr \ M, 



B \ ^y^B 



27r 



Im{Il^+A)) ~ -2ri.T7^ > (35) 



Mb 

It is observed that the twist-3 amplitudes in Eq. ([12|) possess the same power law as the 
r^ terms in Eq. (|lT]). The difference is only the 0(1) ratio M^/A. As stated before, the 
potential linear divergence in the r^ terms alters the naive power counting rules from twist 
expansion. 

Because of 0{Mb/A) ~ tt^ for Mb ~ 5 GeV, the two terms in Eq. ( p4D almost cancel 
each other. The imaginary part in Eq. ( P5| ) is enhanced by the factor 2tt. There is no helicity 
suppression in this case. It is then understood that the annihilation amplitudes from O5 g are 



mainly imaginary, and their magnitude is few times smaller than the factorizable emission 
one as exhibited in Table II. If the mass Mb changes, the relative importance of the real 
and imaginary parts will change. For example, as Mb increases, the cancellation in Eq. ( |3^ 
is not exact, such that the real part becomes larger. On the other hand, the imaginary 
part in Eq. ( P^D decreases. It has been checked that for Mb = 10 GeV, Re{Ffy_^j^-.) is of 
the same order as /m(F.y^^x). We conclude that the smallness of Re{Ffy^^-.) in Table II 
is due to the special value of the B meson mass Mb ~ 5 GeV. For general Mb, it is more 
appropriate to state that annihilation amplitudes from O^^q, without distinguishing their 
real and imaginary parts, scale like 2ri^/(AMB) according to Eq. (plf). 

The above reasoning is applicable to nonfactorizable amplitudes. It can be found, referred 
to Eqs. (^91) -(pip and to the asymptotic behavior of the meson distribution amplitudes, that 
the twist-2 term of each nonfactorizable diagram scales like 1/{AMb). However, because of 
the soft cancellation between a pair of nonfactorizable diagrams in the end-point region of 
X3, the sum of twist-2 terms turns out to scale like l/M"^. The twist-3 and Olr"^) terms in 
each nonfactorizable diagram scale like 

r-K r-KT^ M^ 1 

- ^K^TT^r ' (36) 



AMb ' A2 "' a AMb 

respectively. The cancellation between a pair of nonfactorizable diagrams modifies the above 
power behaviors into 

Ml ' AMb '"'' a M| ' ^ ' 

respectively. For the nonfactorizable annihilation amplitudes (see Eqs. (|53|) and (0)), the 
soft cancellation does not exist, since the B meson is a heavy-light system. However, it 
can be shown that the twist-2, twist-3 and 0{r'^) terms in each nonfactorizable annihilation 
diagram possess the power behaviors, 

J_ Jj^ M^±_ /OCX 

M2 ' M2 ' """^ A M| ' ^^ 

respectively. 
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We emphasize that it is more appropriate to count the power of each individual dia- 
gram, instead of the power of sum of diagrams. In some cases, factorizable contributions 
are suppressed by a vanishing Wilson coefficient, so that nonfactorizable contributions be- 
come dominant. For example, factorizable internal- IV emisson contributions are strongly 
suppressed by the Wilson coefficient 02 in the B -^ J/tpK^*^ decays [Q. In some cases, 
such as the B — > Dn decays, there is no soft cancellation between a pair of nonfactorizable 
diagrams, and nonfactorizable contributions are significant [§. In summary, we derive the 
relative importance of the various topologies of amplitudes given in Eq. (|^). The annihi- 
lation and nonfactorizable amplitudes are indeed negligible compared to the factorizable 
emission ones in the heavy quark limit. However, for the physical mass Mb ~ 5 GeV, the 
annihilation contributions should be included. 

B. Factorization Formulas 

Below we calculate the hard amplitudes for the emission and annihilation topologies, 
whcih have been obtained independently in [|18|. The factorizable penguin contribution is 

written as 

FfW = -SttCfMI / dxidx3 / hidhih^dh^B{xiM) 

Jo JQ 

X I [(1 + X3) ^K (X3) + rx(l - 2x3) ($^ (X3) + $^ (X3))] 

xE(^)(t«)/le(xi,X3,&l,63) 

+2rx$^ (X3) Ei'^ (tf ) h, (X3, xi, 63, bi) } . (39) 

The factorizable annihilation contribution is given by 

F['^' = fS^^ + F.?^^ , (40) 



with 



F^^^'^^ = SttCfMI [ dx2dx3 H h2dh2hdh:i 
Jo Jo 
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x| [(1 -X3)$<^(a;2)$x(x3) 

+2r^r^$; (X2) ((2 - 0:3)$^ (X3) + ^3$^ (X3)) 

- [X2(^^ {X2) $X (X3) 

+2rKr^ ((1 + X2)% (X2) - (1 - X2)$^ (X2)) <l>^ (xs) 
x4J(t(2))/,„(l_a;3,a;2, 63,^2 



(41) 



a6 



+ 



SnCpM'^ / dx2dxs / h2dh2h^dh^ 
Jo Jo 

I [2rx(l - X3)<I'<^(X2) ($^(X3) + $^(X3)) + 4r^% {X2) ^k (xs) 
E'^^{t^^^)ha{x2A-X:,MM) 
4rx$0(x2)$^(x3) + 2a;2r^ {%{x2) - $$, (0:2)) ^k{x^) 

4?Oi'0^a(l-X3,X2,63,&2 



(42) 



for the light quarks q = u and d. Fa in Eqs. (^TJ) and ( P^D are the same as -^^4 but with 
the Wilson coefficient 04 replaced by 02. 

The factors E{t) contain the evolution from the W boson mass to the hard scales t in 
the Wilson coefficients a{t), and from t to the factorization scale 1/b in the Sudakov factors 

s{ty. 



Ei'^^t) =as{t)a'f\t)SB{t)SK{t) 
Ei^\t) =asit)af'\t)S^it)SKit). 



(43) 



The hard functions /I's are 



he{Xi,Xs, 61, 63) = Kq (y/XiX'sMBbi) Stixs) 

X [e{h - b,)Ko (v^Mb6i) Jo (V^sMBbs) 
+e{bs-bi)Ko{^3MBbs)Io{V^3MBbi)] , 

2 
ha{x2,X3,b2,b3) = fy] H^ {V^2X3MBb2) Stixs) 



(44) 
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+^(63 - h)H^o^ iV^sMnh) Jo (v^Msfes)] , (45) 

where St{x) is the evolution function from threshold resummation discussed in Sec. II, and 
Ko, Iq, Hq and Jq are the Bessel functions. 

The hard scales t are chosen as the maxima of the virtualities of internal particles involved 
in the hard amplitudes, including 1/6j: 



t^P = max(V^MB, 1/61, 1/6; 



'3; 



tf = max(v/^MB, 1/h, I/63) , (46) 

t« = maxiVT^TsMB, I/62, 1/&3) , 

t(2) = max(y^MB, 1/62, 1/&3) , (47) 

which decrease higher-order corrections [^. The Sudakov factor in Eq. (|^) suppresses 
long-distance contributions from the large b {i.e., large as{t)) region, and improves the 
applicability of PQCD to B meson decays. We emphasize that the special intermediate scales 



t ~ 0{yAMB) lead to predictions for penguin-dominated decay modes, such as i? ^ (f)K, 
which are larger than those from the factorization and QCD factorization approaches. When 
PQCD analyses are extended to 0{a^) |^, the hard scales can be determined more precisely 
and the scale independence of our predictions will be improved. The 0(0^) corrections to 
two-body nonleptonic B meson decays have been computed in the generalized factorization 
approach [^,^], which indeed improve the scale independence of the predictions. 

For the nonfactorizable amplitudes, the factorization formulas involve the kinematic 
variables of all the three mesons [|^, and the Sudakov factor is given by S" = SbS<s,Sk- 
Their expressions are 

M^ = m:^'^ + m:}"^ + m:^'^ + m:^'^ , (48) 

with 

A^fi^^ = IQ-nCpMlJm, f d[x] r hidhih2dh2^B{xiM)^^{x2) 

* JO JO 
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X| [{X2 + X3)<^k{x3) - rKX3 (^^(xs) + ^^(xs))] 

- [(1 - X2)^k{x3) - rKX3 {^1{X3) - ^^(Xs))] 
X^ef {tf) hf (Xi, X2, X3, 61, &2 



M:t' 



IQ-nCFMl^flN.j d[x] p bidbib2db2<^B {xi, 61) $^ (xa) 

x| [x2^Kix3) - rKX3 (^^(xg) - ^K^^z))] 

xE^£{tf)h^P{x^,X2,X3MM) 

- [(1 - a;2 + X3)^k{.X3) - rKX3 (^^(Xg) + $x(x3))] 

X^ef (4^^) ^?^ (Xl, Xa, X3, 61, 62 



A^2«) = -lQTrCFMlJ2N, f d[x] /°°6irf6iM&2$B (xi, 61) 

* JO Jo 

x| [r^xs ($$,(0:2) - %{X2)) <^k{x3) 

+rKr^ {X2 - X3) ($$, (3:2) ^^ (xs) + $; (X2) $^ (X3) 

+rKr<p{x2 + 0:3) ($$, (3:2) ^K (3^3) - $^ (X2) ^^K (^3) 

X^ef(4^^) tf ^(a^l, ^2, Xg, &1, 62) 



+ 



(49) 



(50) 



r^(l - 0:2) (<I';(X2) + %{X2)) ^k{x3) 

+rKr^{l -X2- X3) ($$, (0:2) $^ (xs) - % [X2) ^K (3^3) 

-rKr4>(\ -X2 + X3) ($$, (X2) $x (3^3) - $0 (X2) $^ (xg) 

xE(f(tf)/.f(a;i,a;2, 0:3,61, 62)}, (51) 

where Nc is number of colors and d[x\ denotes dxidx2dx3. The amphtude M.^^ is the same 
as M-eA but with the Wilson coefficient 0/3 . The nonfactorizable annihilation amplitudes 
are given by 



M^^"^ = M 



Pii) 

a3 



M 



Pii) 

a5 ; 



(52) 



with 
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A<f(^) 



'a3 



-l6nCFMlJ2Nc [ d[x] H hidhih2dh2^B{xiM) 

* JO JQ 

x| [(1 - Xg) $^ (Xz) $i^ (Xg) 

+rKr^ (1 - a;2 - xs) ($$, (^2) ^^ (^3) - $; (x^) ^^ 
-rxr<^(l + xa - 0:3) ($^(a;2)<l'^(a;3) - ^'Axi)^'^^^^^) 



xE^£' (tf^) hf (a;i, a;2, x-^MM) 



- [a;2$^(a;2)$i^(a;3) 

-rxr<^(l -X2- X3) ($^(x2)$^(a;3) + $^(x2)<l'^(a;3) 

+rKr^{l -X2 + X3) (<l'^(x2)<l'^(a;3) + (3 + X2 - x^)%{x2)^'Ki.x^) 

XEif{tf)hf\x,,X2,Xs,h,b2 



Ma^"^ = -IQ-KCpMlxflNc I d\x\ r bidhib2dh2^B{xiM] 

^ Jo Jo 

x| [r^X2 ($*^ (3:2) + % (X2)) ^xixs) 
-rKil - xs)^4x2) (Kixs) - $^(^3))] 
x^if ' (t^^^) /ij:^^ (a;i, 0:2, ^3, 61, 62) 



+ [r^(2 - a;2) (<l>;(a;2) + <l>^(a;2)) <^k{x3) 

-TKil + Xs)<!>^{x2) ($^(X3) - $^(X3))] 

x^if ' (tf) hf (a;i, a;2, 3:3, 61, 62 

The amplitude A^a is the same as Ad^^ but with Wilson coefficient a'^. 
The evolution factors are given by 

E^f (t) =as{t)af{t)S{t)\,,=,,, 
E^ht) =a,(t)ap(t)5(t)|,3=,,. 

The hard functions W\ j = 1 and 2, are written as 

/i? = [e{h - b2)Ko {DMbW) h {DMBh2) 

+e{h2 - h)Ko {DMBh2) /o {DMBhi)] 



(53) 



(54) 



(55) 
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xKo{DjMBb2) , 



x'-^hP 



D]\MBb2 



for D? > 



for D]<0, 



h^'^ 
"/ 



with the variables 






e{h - h2)H^Q~^ (FMBh) Jo {FMBh2) 
+e{b2 - 6i)iJi') (FMBh) Jo {FMBhy 



xKoiFjMBh) , 



in (-1) 

X —^0 
2 ° 



Ff\MBh 



for F/ > 



for F/ < 



(56) 



(57) 



Dl = {xi - a;2)a;3 , 

F>1 = -(1 -xi -a;2)x3 , 

F^ = Xs (1 - Xg) , 

F^ = {xi - X2) (1 - 2:3) , 



F2 = xi + a;2 + (1 - xi - X2) (1 - X3) 



The hard scales t'--'^ are chosen as 



t^P = max (^DMb, \/\Df\MB, l/h, I/62 
tf = max fz^Mfi, vI^Mb, l/h, I/62 



t?) = max (FMb,v/|F2|Mb, 1/61, 1/62) , 



tf) = max ( FMb, \/\Fi\MB, I/61, 1/62 



In the above factorization formulas the Wilson coefficients are defined by 

_ ^ I C'2 / _ ^"1 



.(<?)' 



Nr 



^{^3 + -eqCo , 
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(58) 



(59) 



(60) 



M) 






4" 



05 



"6 






C. 



a 



Ce + ^^ + -eq[Cs + ^r 






Nr 



a(^) = a^^) + ai'^ + af 



IV. NUMERICAL ANALYSIS 



For the B meson distribution amplitude, we adopt the model 
$s(x, h) = Nbx'^{1 - x)^ exp 



1 fxMB\^ ujW 



B'- 



(61) 



2\ ujb J 2 

with the shape parameter ujb = 0.4 GeV. The normalization constant Nb = 91.784 GeV 
is related to the decay constant /b = 190 MeV (in the convention /^ = 130 MeV). As 
stated before, $5 has a sharp peak at x ~ A/Mb ~ 0.1. The and K meson distribution 
amplitudes are derived from QCD sum rules ||36|j50[| , 

3/^ 



$0(x) 

$i (x) 



rT 



2^2N, 



x[l — x) , 

3(1 - 2xf + 0.21 [3 - 30(1 - 2xf + 35(1 - 2xY 



-0.69fl + (l-2x)ln— ^) 
V 1 — xj 



^k{.x) 



^^(x) 



Id) 



4v/2]V; 

v/2iVe 

Ik 



3(1 - 2x)(4.5 - 11.2x + 11.2x2) + 1.38 In 



X 



1 — X 



X 



2^j2Nl 
Ik 

2^/2N^ 



1 - x) {1 + 0.51(1 - 2x) + 0.3[5(1 - 2x)2 - 1]} , 
"1 + 0.24^2/^(1 - 2x) - 0.11Cl/^(l - 2x)l , 
(1 - 2x) [1 + 0.35(10x2 - lOx + 1) 



(62) 

(63) 
(64) 
(65) 
(66) 

(67) 
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with the Gegenbauer polynomials 

To derive the coefficients of the Gegenbauer polynomials, we have assumed Mk = 0.49 
GeV and ttiq = 1-7 GeV. The terms 1 — 2x, rendering the kaon distribution amplitudes 
a bit asymmetric, corresponds to the SU{'i) symmetry breaking effect. We employ Gp = 
1.16639 X 10-^ GeV-2, the Wolfenstein parameters A = 0.2196, A = 0.819, and Rb = 0.38, 
the unitarity angle 03 = 90°, the masses Mb = 5.28 GeV and M<j> = 1.02 GeV, the decay 
constants /^ = 237 MeV, /J = 220 MeV and Jk = 160 MeV, and the 5)] {B-) meson 
lifetime r^o = 1.55 ps (r^- = 1.65 ps) ||5l|. Note that the B —>■ (j)K branching ratios are 
insensitive to the variation of 03. 

We present values of the factorizable and nonfactorizable amplitudes from the emission 
and annihilation topologies in Table III. Contributions from twist-2 and two-parton twist-3 
distribution amplitudes are displayed separately. It is found that the latter, not power- 
suppressed, are in fact more important for f\j,F^. According to the power counting in Sec. Ill, 
the twist-2 contributions to the annihilation amplitudes fsF^ are negligible. This has been 
explicitly confirmed in Table III. As expected, the factorizable amplitudes f^F^ dominate, 
and the annihilation amplitudes fsF^ are almost imaginary and their magnitudes are only 
few times smaller than f,j,F^. The nonfactorizable amplitudes M^ and M^ are down by a 
power of K/Mb ~ 0.1 compared to the factorizable ones f^F^ and /bF^, respectively. The 
cancellation between the twist-2 and twist-3 contributions makes them even smaller. Ma 
and /s-Fa from the operators O12 are of the same order because of the partial cancellation 
between the two terms in the factorization formula for Fa (helicity suppression). 

We demonstrate the importance of penguin enhancement in Table IV. It has been known 
that the RG evolution of the Wilson coefficients 6*4^6 (^) dramatically increases as t < mfy/2, 
while that of Ci,2(t) almost remains constant ||5^. With this penguin enhancement of about 
40%, the branching ratios of the B -^ Kit decays, dominated by penguin contributions, are 
about four times larger than those of the B -^ tttt decays, which are dominated by tree 
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contributions. This is the reason we can explain the observed B -^ Ktc and vrvr branching 
ratios using a smaller unitarity angle 03 < 90° 0,^. In the factorization approach |11[] and 



in the QCD factorization approach ||T2[, it is assumed that factorizable contributions are 
not calculable. The leading contribution to a nonleptonic decay amplitude is then expressed 
as a convolution of a hard part with a form factor and a meson distribution amplitude. In 
both approaches the hard scale is rrih and the intermediate scale AMb can not appear, so 
that the dynamical enhancement of penguin contributions does not exist. To accommodate 
the B — > Kit data in the factorization and QCD factorization approaches, one relies on 
the chiral enhancement by increasing the mass ttiq to a large value rriQ ~ 3 GeV, or on a 
large unitarity angle 03 ~ 120° |T^, which leads to constructive (destructive) interference 



between penguin and emission amplitudes for the B -^ Kn {B -^ titt) decays. 

Whether dynamical enhancement or chiral enhancement is essential for the penguin- 
dominated decay modes can be tested by measuring the B — > (pK modes. In these modes 
penguin contributions dominate, and their branching ratios are almost independent of the 
angle 03. Since is a vector meson, the mass ttiq is replaced by the meson mass M^ ~ 1 
GeV, and chiral enhancement does not exist. Annihilation contributions can not enhance 
the B — i> (f)K branching ratios too much, because they are assumed to be a l/rrih effect in 
the QCD factorization approach [l^. In the PQCD approach annihilation amplitudes reach 



40%, which is reasonable according to our power counting. However, they, being mainly 
imaginary, are not responsible for the large B -^ (pK branching ratios as shown in Table 
IV. If the B -^ (f)K branching ratios are around 4 x 10"^ |ll3| , p!4| , the chiral enhancement 



may be essential. If the branching ratios are around 10 x 10~^, the dynamical enhancement 
may be essential. Therefore, the B -^ (pK decays are the appropriate modes to distinguish 
the QCD and PQCD factorization approaches. The branching ratios of B^ -^ cpK^ and of 
B"^ -^ 0i^^ are almost euqal. We have also evaluated the CP asymmetries of the B -^ (pK 
decays, and found that they are not significant: their maxima, appearing at 03 ~ 90°, are 
less than 2%. 

We emphasize that mo(/i), appearing along with the twist-3 kaon distribution amplitudes, 
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is defined at the factorization scale 1/6 as low as 1 GeV in the PQCD formalism P,^. 
Hence, its value should be located within 1.6 ± 0.2 GeV ||3^. Between the hard scale and 
the factorization scale, there is the Sudakov evolution. In the QCD factorization approach, 
molfi) defined at nib is as large as 3 GeV, leading to chiral enhancement. It has been argued 
that rriQ^jj) and aQ^jj) form a scale-independent product (that is, rriQ^fi) increases, while aQ^fi) 
decreases with /i), such that chiral and dynamical enhancements can not be distinguished 
53| . However, dynamical enhancement exists in both twist-2 and twist-3 contributions, but 



chiral enhancement exists only in twist-3 ones (see Eq. (|39|)). Therefore, they are indeed 
different mechanism. 

At last, we examine the uncertainty from the variation of the hard scales t, which provides 
the informaiton of higher-order corrections to the hard amplitudes. We notice that this is 
the major source of the theoretical uncertainty. The values of ub and rriQ have been fixed at 
around 0.4 GeV and 1.7 GeV, respectively, which are preferred by the B -^ Ktt, vrvr data. 
The light meson distribution amplitudes have been fixed more or less in QCD sum rules. 
The possible 30% variation of the coefficients of the Gegenbauer polynomials lead to minor 
changes of our predictions. Since the analyses of the i? — > tt and B ^ K form factors are 
the same, we constraint the ranges of the hard scales t, such that our predictions for the 
B -^ Ktt branching ratios are within the data uncertainties. The resultant approximate 
range of the hard scales te is given by 

max(0.75v/i^MB, l/h, I/63) < t^ < max(1.25v/i^MB, l/h, 1/63) , 
max(0.75v/^Mij, I/61, 1/63) < tf^ < max{1.25^MB, I/61, 1/63) • (69) 



Note that the coefficients of the factorization scales 1/6, associated with the definition of 
the meson distribution amplitudes, do not change. The variation of the other hard scales 
t is similar, but does not affect the results very much. The theoretical uncertainty for the 
B -^ (j)K branching ratios in Eq. (0) is then obtained. 
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V. CONCLUSION 

In this paper we have shown that a leading-power PQCD formahsm should contain 
contributions from both twist-2 and two-parton twist-3 distribution amplitudes. Threshold 
and kr resummations are essential for infrared finite PQCD analyses of B meson decays. 
Without Sudakov suppression from these resummations, all topologies of decay amplitudes 
possess infrared (logarithmic or linear) divergences. We have explained the power counting 
rules for the factorizable (also nonfactorizable) emission and annihilation amplitudes under 
the sufficiently strong Sudakov effects. The annihilation and nonfactorizable amplitudes are 
suppressed by 2mo/MB and by A/Mb in the heavy quark limit, respectively, compared to the 
factorizable emission ones. For the physical mass Mb ~ 5 GeV, the former should be taken 
into account. In the PQCD formalism the annihilation amplitudes can be calculated in the 
same way as the emission ones without introducing any new free parameters. Hence, our 
formalism has more precise control on the annihilation effects than the QCD factorization 
approach. Annihilation contributions of 40% at the amplitude level are reasonable according 
to our power counting. However, these amplitudes are not responsible for the large B -^ (pK 
branching ratios, since they are mainly imaginary. 

We have emphasized that exclusive heavy meson decays are characterized by a lower 
scale AMb, for which penguin contributions are dynamically enhanced. This enhancement 
renders penguin-dominated decay modes acquire branching ratios larger than those from the 
factorization and QCD factorization approaches, even when the final-state particle is a vector 
meson. We have proposed the B — > (j)K decays as the ideal modes to test the importance 
of this mechanism. If their branching ratios are as large as 10 x 10^® (independent of the 
unitarity angle (ps), dynamical enhancement will gain a convincing support. The answer will 
become clear, when the consistency among the BaBar, Belle, and CLEO data is achieved. 
We have also found that the CP asymmetries in the B -^ (pK modes are vanishingly small 
(less than 2%). 
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APPENDIX A: TWO-PARTON DISTRIBUTION AMPLITUDES 

1. B Meson Distribution Amplitudes 



The B meson distribution amplitudes are written as [1C,54] 



with the dimensioless vectors n+ = (1,0, Oj-) and n_ = (0,1, 0^)- As shown in |10], the 
contribution from $b is negligible, after taking into account the equation of motion between 
$B and $B- Hence, we consider only a single B meson distribution amplitude in the heavy 
quark limit in this work. As the transverse extent b approaches zero, the B meson distribu- 
tion amplitude ^b{x, h) reduces to the standard parton model ^b{,x) = ^b{x, 6 = 0), which 
satisfies the normalization 

2. (j) Meson Distribution Amplitudes 

To define the meson distribution amplitudes, we consider the following nonlocal matrix 
elements 
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miA^mP2)) = uM^ 



p, 



2/x 



e ■ z 
'K~zJo 



t/xae-"^^^" 



"011 (a;2) 



Jo 



X2 



e ■ z 



{P2-Z) 



-Mi 



dx2e-'^''^^-' 



{0\m^,MzmP2))=^f^ 



[^TfiP2u — ^TuP'. 



+ {P2iJ.Z,y — P2iyZfj 



2fi) 



e ■ Z 



1 

2' 



(P2 ■ ^)2 



g3[X2) 



dx2e-''''^^-\ 

Ml / dx2e 




(A3) 



'Tia;2J 



-To ■ 2; JO 



' / «T^ « ^(iLci 



\rsiO)IsizMiP2)) = -[fl~U 



zMi 



dx2e-'"^^^"/i3(x2) 



cixae-^^^^-^/ii'Vxsl 



(A4) 
(A5) 



2 V^ ^'"M,^/ wo 

where /^ and /J are the decay constants of the (p meson with longitudinal and transverse 
polarizations, respectively, er the transverse polarization vector, X2 the momentum asso- 
ciated with the s quark at the coordinate z ex (0,1, 0^), and mg the s quark mass. The 
explicit expressions of the distribution amplitudes 0, g and h with unity normalization are 
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The contributions from the distribution amplitudes gj! , (J)t and h^ vanish for two-body 
B meson decays, in which only longitudinally polarized mesons are produced. The con- 
tributions from 0||, h» , h\f and (?3 are twist-2, twist-3, twist-3, and twist-4, respectively. It 
is easy to confirm that g^ does not contribute to the factorizable emission and annihilation 
amplitudes. The term proportional to the small ratio 2ms /M^ in Eq- ( [A5| ) is negligible. 
Therefore, up to twist 3, we consider the following three meson final-state distribution 
amplitudes. 






$</,(a;2) 



^ /^2 -.t 



V2N, 



ni^2) 



Ma 



V2N, 



■m^2) 



(A6) 



with 
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Jd 



-.h 



it) 



$! 



fl d 



M 



Ay/2N^dx 



(A7) 



2v/2i\r: " "^ 2v/2Ar: 

The spin structures associated with the meson distribution amplitudes can be derived 
from Eqs. ([A3D-( 
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3. Kaon Distribution Amplitudes 

The general expressions of the relevant nonlocal matrix elements for a kaon are given by 

Jo 
-^kMI-^ f dx,e-'-^''-^gK{xs) , (A8) 

{^\mi^n{z)\K{P^)) = -zfKmo t dx^e-'^^'^'-'M^s) , (A9) 

Jo 



i „ f. Ml' 



(0|5(0)75a^.M(z)|i^(P3)) = 7.//^mo U - — r {Pz^.z, - P-^.z^) 



6 y rriQ ^ 

X f' dx3e-'^^''^-'Mx3) , (AlO) 

Jo 



with the mass 



mo = -^^^ . (All) 






fx is the kaon decay constant and x^ the momentum fraction associated with the u quark 
at the coordinate z oc (1, 0, Oy). The explicit expression of the wave fucntions and qk with 
unit normalization are referred to p6| . 

The contributions from the distribution amplitudes 0„, 0p, (j)fj and gx are twist-2, twist-3, 
twist-3, and twist-4, respectively. Note that g^ does not contribute to factorizable emission 
and annihilation amplitudes. Hence, the factorizable annihilation amplitude in Eq. ( ^T] ) 
is complete in the r^ terms. We consider the following three kaon final-state distribution 
amplitudes 

with 

where the term Ml/rn^ in Eq. ( |A10D has been neglected. The spin structures associated 
with the kaon distribution amplitudes can be derived from Eqs. (|A^)- (|A10|) . 
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Figure Captions 

FIG. 1: Leading contribution to the B -^ (pK decays, where H denotes the hard amplitude 
and ki, i = 1, 2, 3 are the parton momenta. 

FIG. 2: Lowest-order diagrams for the B ^ ir,K transition form factors. 



37 



TABLES 

TABLE L B ^ iT,K transition form factors without threshold resummation (column A) and 
with threshold resummation (column B). 



Form factors 


F^^(O) 


F^^(O) 


uJB (GeV) 


A 


B 


A 


B 


0.35 


0.603 


0.356 


0.703 


0.430 


0.36 


0.579 


0.343 


0.675 


0.413 


0.37 


0.557 


0.330 


0.647 


0.398 


0.38 


0.535 


0.318 


0.621 


0.382 


0.39 


0.515 


0.306 


0.597 


0.368 


0.40 


0.496 


0.295 


0.574 


0.354 


0.41 


0.478 


0.285 


0.552 


0.342 


0.42 


0.461 


0.275 


0.532 


0.329 


0.43 


0.445 


0.266 


0.512 


0.318 


0.44 


0.429 


0.257 


0.494 


0.307 


0.45 


0.414 


0.248 


0.476 


0.296 
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TABLE II. Time-like form factors F'f'^ 



Form factors 


F^^(O) 


F^^(O) 


^(t% xlO^ 


p<t>^ X 1 02 


mo(GeV) 


Re 


Re 


Re 


Im 


Re 


Im 


1.4 


0.295 


0.312 


1.49 


0.56 


-2.58 


12.72 


1.5 


0.308 


0.326 


1.61 


0.57 


-3.06 


13.13 


1.6 


0.320 


0.340 


1.71 


0.62 


-3.20 


13.41 


1.7 


0.333 


0.354 


1.78 


0.63 


-3.48 


13.52 


1.8 


0.345 


0.368 


1.92 


0.68 


-3.79 


13.88 


1.9 


0.357 


0.383 


2.02 


0.73 


-3.99 


14.67 


2.0 


0.370 


0.396 


2.10 


0.74 


-4.40 


14.82 
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TABLE III. Twist-2 and higher-twist contributions to the B -^ (j)K decay amphtudes. 



Decay 


B^ -^ (t)K^ 


Components 


twist-2 (] 


LO-4 GeV) 


higher-twist (10""^ GeV) 


total (10 


-4 GeV) 


Amphtudes 


Re 


Im 


Re 


Im 


Re 


Im 


f<t>Fe 


9.48 


— 


27.71 


— 


37.19 


— 


Mf 


-3.30 


2.44 


1.60 


-1.11 


-1.70 


1.33 


IbFP 


0.07 


-0.01 


-2.57 


-15.47 


-2.50 


-15.48 


Mf 


0.06 


0.41 


0.30 


-0.25 


0.36 


0.16 


fsFa 


-1.21 


0.47 


41.00 


13.52 


39.79 


13.99 


Ma 


0.88 


-12.74 


-7.05 


2.61 


-6.17 


-10.13 



Decay 


B^ ^ 4>K'' 


Components 


twist-2 ( 


10-^ GeV) 


higher-twist (10-^ GeV) 


total (10 


-^ GeV) 


Amplitudes 


Re 


Im 


Re 


Im 


Re 


Im 


UF^ 


9.48 


— 


27.71 


— 


37.19 


— 


MP 


-3.30 


2.44 


1.60 


-1.11 


-1.70 


1.33 


IbFP 


0.07 


-0.01 


-2.68 


-15.76 


-2.61 


-15.77 


Mf 


-0.03 


0.77 


0.48 


-0.28 


0.45 


0.49 
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TABLE IV. Enhancement effects in the B^ -^ (j)K^ decay amphtudes 



Scales 


^ = t 


fi = 2.5 GeV 


Amplitudes 


Re (10-4 GeV) 


Im (10-4 GeV) 


Re (10-4 GeV) 


Im (10-4 GeV) 


UFf 


37.19 


— 


23.14 


— 


Mf 


-1.70 


1.33 


-1.05 


0.62 


fsFr 


-2.50 


-15.48 


1.92 


-12.83 


M^ 


0.36 


0.16 


-0.05 


0.19 


fsFa 


39.79 


13.99 


37.73 


13.14 


Ma 


-6.17 


-10.13 


-0.56 


-9.05 


Branching Ratio 
(without Ann.) 


9.8 X 10-6 


3.8 X 10-6 


Branching Ratio 
(with Ann.) 


10.2 X 10-6 


5.6 X 10-6 
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